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We study photon correlations generated by scattering from three-level systems (3LS) in one dimension. The 
two systems studied are a 3LS in a semi-infinite waveguide (3LS plus a mirror) and two 3LS in an infinite 
waveguide (double 3LS). Our two-photon scattering approach naturally connects photon correlation effects with 
inelastically scattered photons; it corresponds to input-output theory in the weak-probe limit. At the resonance 
where electromagnetically induced transparency (BIT) occurs, we find that no photons are scattered inelastically 
and hence there are no induced correlations. Slightly away from BIT, the total inelastically scattered flux is large, 
being substantially enhanced due to the additional interference paths. This enhancement carries over to the two- 
photon correlation function, which exhibits non-classical behavior such as strong bunching with a very long 
time-scale. The long time scale originates from the slow-light effect associated with BIT. 


I. INTRODUCTION 

The many similarities between quantum transport of elec¬ 
trons (conduction) and optical phenomenon (propagation of 
EM radiation) have been used over the years to enrich both 
fields. While a scattering approach to the propagation of 
light, with input and output amplitudes, is quite natural in 
both classical and quantum optics [1, 2], a comparable ap¬ 
proach to electronic phenomena developed slowly. First in¬ 
troduced by Landauer [3, 4], it was subsequently substantially 
developed by Biittiker [5-7]. This approach was then used, 
for instance, to develop parallels in mesoscopic physics be¬ 
tween electronic and photonic phenomena, such as coherent 
backscattering of electrons or photons from disordered me¬ 
dia [8, 9]. Another example is in the development of semi- 
classical (or eikonal) approximations to quantum chaotic phe¬ 
nomena and the inclusion of diffractive effects [10]. While 
these parallels were developed mainly in the non-interacting- 
particle or linear-optics regime, interacting particles and the 
corresponding nonlinear regime are, of course, of key interest 
in both photonic and electronic transport. One particular set¬ 
ting that has received a great deal of attention in the quantum 
transport community is one-dimensional (ID) electrons inter¬ 
acting with local quantum impurities, a setting that includes 
for instance the Kondo problem, Anderson impurity model, 
and Bethe-Ansatz solutions [11-13]. The parallel photonic 
system is a one-dimensional EM waveguide strongly coupled 
to discrete non-linear quantum elements such as atoms, quan¬ 
tum dots, or qubits; in analogy with “cavity QED” [2, 14], the 
study of such systems is known as “waveguide QED.” 

The study of waveguide QED has increased rapidly over the 
past decade. Prior to that, there were a few early papers on the 
subject [15-19] that, for instance, exploited many-body ap¬ 
proaches developed for electronic problems. The dramatic in¬ 
crease in interest starting in the period 2005-2008 [20-24] was 
driven by experimental progress toward achieving strong cou¬ 
pling between the waveguide and the local quantum system. 
Indeed, several experimental waveguide-QED platforms are 
being actively pursued. These include superconducting qubits 
coupled to a microwave transmission line [25-30], semicon¬ 
ductor quantum dots coupled to either a metallic nanostructure 


[31, 32] or a photonic-crystal waveguide [33], and more tra¬ 
ditional quantum optics settings in which atoms provide the 
local quantum system and the waveguide is an optical fiber or 
glass capillary [34, 35]. Interesting waveguide-QED effects 
occur when the coupling to the waveguide dominates other 
emission or dephasing processes. Experiments in this inter¬ 
esting regime have been performed in several of the above 
waveguide-QED platforms. 

Two aspects of waveguide QED have attracted particular at¬ 
tention theoretically; the manipulation of single photons and 
the production of non-classical light. In the single photon 
arena, a variety of devices have been proposed that build on 
the manipulation of single photons by qubits or three-level 
systems (3LS) that is possible in ID systems; for represen¬ 
tative work in this area see Refs. [33, 34, 36] and references 
therein. With regard to non-classical light, the main character¬ 
istics studied are the photon-photon correlation function (also 
called the second-order coherence [37]) and the photon statis¬ 
tics. The majority of work on these topics has treated a single 
quantum system coupled to the waveguide, where the single 
quantum system is modeled as a two-level system (2LS) or 
the only slightly more complicated driven 3LS (for very re¬ 
cent work along these lines see, for example. Refs. [38-41]). 
Correlation effects in a multi-qubit waveguide have been stud¬ 
ied in a number of recent papers using a variety of techniques 
[42-55]. In most of these, the Markovian approximation is re¬ 
quired in order to simplify the interactions between the qubits 
via the waveguide [42^8]. There are, however, a few non- 
Markovian results [49-55] which have been used to delineate 
the range of validity of the Markov approximation. 

Here we extend our recent results on multiple 2LS waveg¬ 
uide QED [49-51] to the case in which driven 3LS are 
used. We calculate the two-photon wavefunction and focus 
on photon-photon correlations. We find that these correla¬ 
tions are substantially enhanced in systems containing 3LS, 
making them better experimental candidates for further study 
of the non-classical light produced. Furthermore, we find that 
the complexity of the structure enhances the photon-photon 
correlations—they are enhanced by adding additional nonlin¬ 
ear elements (qubits) as well as by simply adding a mirror. 
The photons can be either bunched or anti-bunched depend¬ 
ing on the situation, and we find cases of both strong bunching 
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FIG. 1. Schematic of waveguide QED. (a) Single 3LS coupled to 
a semi-infinite waveguide with qubit-mirror separation a, cde = tuo 
and uja = ljJq — A; (b) Two identical 3LS, separated by distance L, 
coupled to an infinite waveguide. 


and anti-bunching. 

The paper is organized as follows. In the next section, we 
first recap the standard model of waveguide QED and a 3LS 
and summarize our approach to finding the two-photon wave- 
function. Then we present the physical quantities that are cal¬ 
culated, emphasizing the total inelastic scattering as a measure 
of the correlated part of the wavefunction. Results for a single 
3LS are presented in Section 3 as a basis for comparison to the 
more complex structures studied later. In Section 4 we add a 
mirror to the system, thus studying a single 3LS in a semi¬ 
infinite waveguide. Section 5 covers results for two 3LS in 
an infinite waveguide. In the results of both Sections 4 and 5, 
inelastic scattering is enhanced, suggesting more visible cor¬ 
relation effects. Finally, in Section 6 we discuss implications 
of the results and conclude. 


II. MODEL AND OBSERVABLES 
A. Waveguide QED model with multiple three-level systems 

The standard Hamiltonian of waveguide QED [20, 22] con¬ 
sists of a one-dimensional bosonic field that can travel to the 
left or right coupled to N local quantum systems, often called 
simply qubits. For a schematic see Fig. 1. Within the rotating- 
wave approximation, the Hamiltonian in real space reads (tak¬ 
ing h = c = 1) 


H = Hqs — i I dx 


al{x)^aR{x) - al{x)-^ai^{x) 


N 


— oo 
oo 


^ 14 / dx5{x-Xi) aj,(a;)crW-FcrWa„(a;) 

i^l a^L,R 


( 1 ) 


where Ueg = \s)i{g\ denotes the Pauli raising operator of the 
/-th qubit with position Xi and coupling strength Vi, and or l 
denote the annihilation operators of right- or left- going pho¬ 
tons. The corresponding decay rate of the /-th qubit to the 
waveguide is = 214^. Throughout this paper, the coupling 
of all of the qubits is the same, V. In order to assess the maxi¬ 
mum possible non-classical light effects that could be present, 
we focus on the lossless limit. 

The local quantum systems that we consider here are iden¬ 
tical 3LS, Hqs, — ^ 3 LS- Hamiltonian for a A-type 

3LS is 

-ffsLS = Wo|e)(e| +a;^|s)(s| + ^ (|e)(s| -F |s)(e|), (2) 

in which H is the Rabi frequency of the classical driving and 
ujs = We — A with A being the detuning between the driv¬ 
ing frequency and the frequency of the |s) to |e) transition. 
(The frequency corresponding to the ground state is taken to 
be zero.) Finally, we note that a mirror can be introduced as a 
boundary condition when solving for the single-photon wave- 
function [51, 56-58]. 

To construct the two-photon scattering wavefunction, we 
use the Lippmann-Schwinger equation [49-51, 59], in which 
the Pauli raising and lowering operators, deg and dge, are re¬ 
placed by bosonic creation and annihilation operators, and 
b. (A similar approach has been used in the case of two- 
electron scattering [60, 61].) To satisfy the level statistics, it is 
necessary to introduce an additional on-site repulsion U to be 
taken as infinite at the end. For a 2LS, it is known that this ap¬ 
proach correctly gives all measurable quantities [49, 51]. For 
a 3LS, in addition to repulsion for each upper level, an extra 
term has to be added so that the double occupancy can be fully 
ruled out; the repulsion operator V is 

V=^ {blblbebe + blblbsbs + 2b%blbs) . (3) 

Note that the coefficient of the last term is chosen for con¬ 
venience; any coefficient would be canceled out after taking 
U —>■ 00 . Once a proper on-site interaction V is introduced 
for each qubit, the calculation of the two-photon wavefunction 
\ 1 lj 2 ) can be done straightforwardly [50, 51]. 

The two-photon scattering wavefunctions that emerge have 
a common structure (see, for instance. Ref. [62] and refer¬ 
ences therein); an incoming multi-particle plane wave and 
two contributions to the outgoing part— (i) an outgoing multi- 
particle plane wave in which the momenta of the photons are 
just rearranged and (ii) a part that involves a continuum of 
momenta that cannot be written as a simple plane wave. The 
first part is elastic scattering, and the second is inelastic (for 
each individual particle). It is this latter part of the wavefunc¬ 
tion that encodes the correlation between the two photons; it 
has been called the “bound state part” [21, 59, 63-65] because 
as the distance between the particles grows, the wavefunction 
decays exponentially. In scattering theory, it corresponds to 
the two-particle irreducible T-matrix [62]. 

Having defined the model, we now turn to the observables 
that we study. 
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B. Power spectrum (resonance fluorescence), S{uj) 

The power spectrum is defined as the Fourier transform of 
the amplitude correlation function (also called the first-order 
coherence), 

Sa{uj) = J dt {ilj2\a'l{xo)aa{xo + t)\tjj2), (4) 

for a = R or L, corresponding to right-going or left-going 
photons. Sa{uj) gives simply the intensity of outgoing pho¬ 
tons as a function of their frequency. The elastic scattering 
gives a 5-function contribution which we drop; thus, in this 
paper Sa{uj) refers to the inelastic power spectrum. 

C. Total inelastically scattered component, F{kin) 

Since in the two-photon wavefunction, the inelastically 
scattered component is necessarily the correlated “bound 
state” part of the wavefunction and vice versa, the total in¬ 
elastically scattered power provides a measure of the overall 
strength of correlations in the situation studied. This, then, can 
be used as a figure of merit to compare different systems or 
optimize parameters in order to find the largest non-classical- 
light effects. Thus, we define a function of the incoming mo¬ 
menta, 

F{h^,\,hn,2) = J duj [Sr{uj) + Si^iuj)]. (5) 

In our study, we take the two incoming photons to have the 
same energy, and so we consider a function of one variable, 
F(fcin). We shall refer to the maximum value of F{k), de¬ 
noted Fpeak, and the value of k at which this maximum is 
reached, fcpeak- 

D. Photon-photon correlation (second-order coherence), g 2 (t) 

A key signature of non-classical light is that the photons can 
be bunched or anti-bunched in space or time. This is revealed 
through the photon-photon correlation function, g 2 {t) [37]. 
Physically, p 2 {t) is proportional to the joint detection prob¬ 
ability of two photons far from the scattering region that are 
separated by distance ct. It is customarily normalized by tak¬ 
ing the ratio of the joint detection probability to the probability 
of measuring single photons independently; thus, g 2 {t) = 1 
signifies no correlation. As the photons should be uncorre¬ 
lated at very large separations in the systems we study, we 
have the property g 2 {t —)■ c») = 1 [37]. 

For a weak coherent state, the correlation comes from the 
two-particle wavefunction while the normalization is given 
by the single-photon wavefunction (see, e.g.. Refs. [38, 50]), 
yielding 

_ {'>p2\al,[xo)a'^{xQ + t)aa{xo + t)aa{xo)\'tp2) 
\{iJi\ai{xo)aa{xo)\tl^i)\'^ 

We now turn to presenting results for these observables in 
three cases. 



FIG. 2. Single 3LS coupled to an infinite waveguide: single¬ 
photon transmission probability T, total inelastically scattered flux 
F, normalized power spectrum S, and correlation functions g 2 of 
the transmitted (fourth row) and reflected (fifth row) intensity. The 
classical Rabi frequency for the first column is G = r/4 and for the 
second column is F. T and F are functions of the incident frequency 
k, while the power spectrum S and the correlation function g 2 are 
calculated by fixing k = fcpeak, the incident frequency giving rise to 
maximum F. Note that for a single 3LS, Sr — Sl = S. A hori¬ 
zontal line indicating the uncorrelated value 32 = 1 is plotted. The 
transmittance at fe = fcpeak for G/F = 1/4 is r(fcpcak) = 39.3% and 
for fl/F = 1 is 1.8%. (Parameters used: loq = Us = lOOF and 
A = 0.) 


III. SINGLE 3LS 

We start from a single 3LS coupled to an infinite waveg¬ 
uide, a system previously investigated in Refs. [38, 66-74]. 
For two identical photons injected from the left, the results 
are shown in Fig. 2. When each photon is resonant with the 
3LS, k = uJs = ujg — A, electromagnetically induced trans¬ 
parency (FIT) occurs—an interference effect between differ¬ 
ent processes within the 3LS that leads to decoupling of the 
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3LS from the waveguide [75]. The width of the resulting peak 
in transmission is ^ [72, 75]. BIT is clearly seen in 

the single-photon transmitted intensity as a function of input 
wavevector k in panels (a) and (f) of Fig. 2. 

As an illustrative case, we give the expressions for S{uj) 
and F{k) in the Appendix, and plot them in panels (b), (c), 
(g), and (h). F{k) has a zero at the BIT peak (when k = Wg) 
and two sharp peaks nearby. The fact that F{kin = oJs) = 0 at 
perfect transparency is, in fact, true when any number of iden¬ 
tical 3BS are coupled to the waveguide. This means that at this 
resonance, all of the photons are scattered elastically, and that, 
therefore, there is no “bound state” part of the wavefunction. 
This is consistent with the fact that there are no correlations in 
these cases—the correlation function of the transmitted pho¬ 
tons is < 72 (i) = 1 [38, 71]. This phenomenon has been called 
“fluorescence quenching” [76, 77]. 

The surprisingly large peak value of F in panel (b) sug¬ 
gests that there will be large correlation effects because of the 
necessary connection between inelastic processes and photon- 
photon correlation. For comparison, Fpeak for a single 2BS 
is S/ttB: this can be derived either by setting il to zero in 
Bq. (A2) or by multiplying Bq. (D4) of Ref. [5 1 ] by Ftt. Thus, 
when n/r = 1/4, fpeak is about five times larger in the 3BS 
case than for a 2BS. We infer that a 3BS is much more effec¬ 
tive than a 2BS in creating photon correlations. 

To explore the effect of the correlated part of the wavefunc¬ 
tion, we consider the most favorable case; photons injected at 
the energy at which F is maximized, k = fcpeak- In the BIT 
regime (weak driving, fl/B <C 1 ), panel (c) shows that the 
power spectrum has a sharp peak centered at fcpeak, and panels 
(d) and (e) show that the time scale of the correlation function 
is long, ^ 40/r. The slow decay time is associated with the 
time delay r [51, 75], which is itself associated with the in¬ 
verse of the width of the BIT peak. At k = uig, the time delay 
is simply r = 2 r/n^. 

An alternative route to obtaining a large and strong 
correlations with a long decay time is to use multiple 2BS. 
From our previous work [51], we note that in order to obtain 
the large value of Fpeak and slow decay obtained here for a sin¬ 
gle 3BS, one has to use of order ten 2BS—a considerably more 
complicated structure (but one that does not require classical 
driving). 

On the other hand, features typical for a single qubit, re¬ 
gardless of its level structure, coupled to an infinite waveg¬ 
uide are present. For example, the reflected p 2 (t = 0) is al¬ 
ways zero because a single atom cannot reflect two photons at 
the same time. Moreover, the transmitted and reflected power 
spectra are identical [for a single 3BS, see Bq. (Al)]. 

The magnitude of the held driving the 3BS has a large ef¬ 
fect on the correlation properties—compare the first and sec¬ 
ond columns of Fig. 2. First, the magnitude of F clearly de¬ 
creases, implying less correlation. In fact, for 12 7 ^ 0, i^peak 
is a monotonically decreasing function of 12 [see Bq. (A2)]. 
It can be very large in the BIT regime (fl/B ^ 1) and then 
decreases to IG/ttB as f2/r approaches one; curiously, Fpeak 
at 12 = 0—that is, in the 2BS case—is half this large 12 value. 
Second, the spectrum of outgoing photons is broadened into 
a Borentzian form. Finally, the timescale for decay of g 2 (t) 


is considerably shortened. The nature of the correlations at 
12/r ~ 1 is such that there are few special features. Thus, 
using a 3BS with strong driving provides a route to relatively 
simple, less structured correlation effects, a feature that may 
be desirable in some situations. 

We close this section by commenting on the behavior of 
F{k) at large fc. Upon inspecting the expression for F{k) in 
Bq. (A2), we And that it decays as 1/fc"^ for incident frequency 
fc far away from ujg. This form of decay is also true for a single 
2BS (12 = 0). 


IV. 3LS PLUS MIRROR 

To construct the second system that we study, we add a mir¬ 
ror to the right of the waveguide so that photons can enter or 
leave only from the left—see Fig. 1(a) for a schematic—thus, 
the system consists of a single 3LS coupled to a semi-infinite 
waveguide [78]. As a result, the amplitude of the single¬ 
photon reflection is unity, |r(fc)| = 1. The calculation of the 
single- and two- photon eigenstates is explained in detail in 
Ref [51] in the case of a 2LS; the 3LS case is, of course, very 
similar. Here we suppose the 3LS is placed close to the end of 
the waveguide, fcga < 27r where fco = ujq/c = 27r/Ao is the 
wavevector associated with the qubit frequency ujq. We use, 
therefore, the Markovian approximation. 

Results are shown in Fig. 3 for two values of the separa¬ 
tion between the 3LS and the mirror: fcoo = 7r/2 and 7r/4 
which correspond to a = Ao/4 and Ao/8. As for a single 
2LS coupled to a semi-infinite waveguide (see, for example. 
Refs. [30, 51, 57, 58, 79] and references therein), interference 
effects associated with the qubit-mirror separation a have a 
large effect on all measurable quantities. These effects can 
be connected to the complicated response of the poles of the 
reflection amplitude to the placement of the 3LS. For a dis¬ 
cussion of the role of poles in the 2LS case see Ref. [51]; the 
present 3LS case is similar. 

We highlight several features of the results in Fig. 3: 

1. The magnitude of inelastic scattering increases substan¬ 
tially upon adding a mirror: compare the first row of 
Fig. 3 to the second row of Fig. 2. In the BIT regime, 
Fpeak is roughly an order of magnitude larger when a 
mirror is present. Note in addition that for both values 
of o, F = 0 when the injected photons are on resonance 
with the 3LS. 

2. F{k) is symmetric for koa = 'k 12 but is highly asym¬ 
metric for fco a = 7r/4. The latter is connected to the 
asymmetric pole structure. 

3. S[lij) is narrow and largely structureless, centered at 
and symmetric with respect to fcin (see the second row 
in Fig. 3). 

4. Strong initial bunching is found (third row), g 2 {t = 
0) > 5, in contrast to the initial anti-bunching found for 
reflection in an infinite waveguide—^because a photon 
reflected by the mirror can stimulate emission from the 
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FIG. 3. Single 3LS coupled to a semi-infinite waveguide [cf. Fig. 1(a)]: total inelastically scattered flux F, normalized power spectrum S, 
and reflected correlation function g 2 . Qubit-mirror separation a and classical Rabi frequency Q, used for each column are labeled on the top. 
F is a function of the incident frequency k, while the power spectrum S and the correlation function 32 are calculated by fixing k = fcpeak, 
the incident frequency giving rise to maximum F. A horizontal line indicating the uncorrelated value 02 = 1 is plotted. (Parameters used: 
= a;, = lOOr and A = 0.) 


3LS, the reflected g2{0) need not be zero. The bunch¬ 
ing is followed by anti-bunching characterized by a long 
time-scale. 

5. Increasing the classical driving held Q causes a shift in 
fcpeak away from the resonant frequency—which is nat¬ 
ural as the EIT feature becomes broader—as well as a 
sharp decrease in the magnitude F(fcpeak)- The mag¬ 
nitude of the secondary peak in the case kQa = 7r/4 is 
less sensitive to Q than the primary peak. The long time 
anti-bunching is cut-off for increasing Q. 

One can derive the time delay r for a semi-infinite waveg¬ 
uide by differentiating the phase (f>k of the reflection ampli¬ 
tude r{k) with respect to k [51, 80]. Surprisingly, on reso¬ 
nance (k = uig), the time delay has a very simple expression 
T = 2 X 2f /Q^, where f = T [1 — cos(2fcoa)] is the effec¬ 
tive decay rate of the 2LS in front of the mirror [30, 51, 57]. 
This suggests that this system can be seen, roughly speaking, 
as a double 3LS (a real one and its mirror image) coupled to 
an infinite waveguide, a system to which we turn in the next 
section. 


V. TWO 3LS 

The last example that we consider is that of two 3LS cou¬ 
pled to an infinite waveguide. A schematic is shown in 
Fig. 1(b), and the results are in Fig. 4. 


Many of the properties of a two-3LS system are similar to 
either those of a single 3LS or those of a 3LS plus a mirror. 
In the first row of Fig. 4, we recognize, for instance, the FIT 
peak in T(fc). F{k) goes to zero at the FIT peak—at perfect 
transmission there are no correlation effects. We And that the 
widths of the FIT peaks in both cases are smaller than that 
of a single 3LS in an infinite waveguide, which qualitatively 
agrees with the previous studies [72, 75]. (Note that in our 
study the inter-3LS scattering is included.) T and F are sym¬ 
metric for koL = tt/2 (L = Ao/4) and distinctly asymmetric 
for koL = irjA (L = Xq/S). The peak value of F is large, 
comparable to that for a 3LS plus mirror. The fact that a more 
complicated structure, providing more interference paths, pro¬ 
duces a larger T"peak, as in these two examples, suggests that 
it may be a general rule. 

The power spectrum of inelastically scattered radiation 
upon excitation at fcpeak is slightly broader than in the other 
two cases and appears to be formed from overlapping peaks. 
An important difference with respect to the single 3LS case is 
that S{uj) is not the same for the reflected and transmitted ra¬ 
diation, though they are of the same order of magnitude. This, 
then, is the same as was seen for 2LS systems [28, 47, 51]. In 
the 2LS case where a photonic band gap builds up as the num¬ 
ber of 2LS increases, (reflection) is always larger than Sr 
(transmission) around the 2LS transition frequency. However, 
Fig. 4(r) shows this is not be true for multiple 3LS. 

Turning to the correlation functions (rows 5 and 6), we 
see that both transmitted and reflected photons are bunched 
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FIG. 4. Two 3LS coupled to an infinite waveguide [cf. Fig. 1(b)]: single-photon transmission T, total inelastically scattered flux F, nor¬ 
malized power spectrum S, and the correlation functions g 2 for the transmitted (fourth row) and reflected (fifth row) intensity. Quhit-qubit 
separation L and classical Rahi frequency Q used for each column are labeled on the top. T and F are functions of the incident frequency 
k, while the power spectrum S and the correlation function g 2 are calculated hy fixing k = fcpcak, the incident frequency giving rise to max¬ 
imum F. The transmitted, reflected, and total power spectra are shown in red dotted, hlue dashed, and black solid curves, respectively. A 
horizontal line indicating the uncorrelated value ^2 = 1 is plotted. The transmittance at A: = Apeak for the four columns are (from left to right) 
37%, 16.9%, 67.4%, and 54.9%. (Parameters used: ujq = ujg = lOOF and A = 0.) 


at f = 0 in all cases shown. There is clearly a long time-scale 
present, in addition to a short one of order l/T. The long 
time scale is naturally explained in terms of the time delay. 
Upon inspecting the single-photon transmission coefficient of 
N identical 3LS (see for example Ref. [72]), one finds that the 
time delay caused by N 3LS is simply N times the single 3LS 
time delay, r = 27VrTherefore, by increasing the num¬ 
ber of 3LS, the time scale of 52 is proportionally lengthened. 


The effect of increasing the classical driving is very simi¬ 
lar to that for a single 3LS; the peaks in F become smaller and 
shift in frequency, and the correlation functions 52 become 
more normal with strong bunching or anti-bunching followed 
by a decay to 52 = 1- 


VI. CONCLUSIONS 


We have studied non-classical light produced in two 3LS 
systems—a 3LS in a semi-infinite waveguide and two 3LS in 
an infinite waveguide—by extending the formalism that was 
previously applied to multiple 2LS. Sharp peaks in the inelas¬ 
tic flux and the power spectrum appear near the two-photon 
resonance at which EIT exists. By tuning the classical driving 
or the qubit-mirror or qubit-qubit separation, these peaks can 
be made narrower and higher, indicating stronger correlations 
or “bound state” effects. The photon-photon correlation func¬ 
tion is studied at the frequency at which the inelastic flux is 
maximized; it shows complicated bunching and anti-bunching 
with a very long time-scale. Our study reveals how the slow- 
light effect associated with EIT is expressed in higher-order 
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quantities such as the correlation function. 

Our results show that a 3LS is much more effective than 
a 2LS in creating photon correlations, making them better 
candidates for experimental study of the non-classical light 
produced. Furthermore, we find that the complexity of the 
structure enhances the photon-photon correlations—they are 
enhanced by adding additional nonlinear elements (qubits) as 
well as by simply adding a mirror. 
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Appendix A: Power spectrum of single 3LS without a mirror 

To check the validity of our scattering approach [49-51], 
let us focus on a single 3LS coupled to an infinite waveg¬ 
uide, since it is widely known that the wavefunction of pho¬ 
tons propagating in a ID waveguide coupled to a single, lo¬ 
cal quantum system can be solved exactly without much dif¬ 
ficulty. For this case, we compared our two-photon scattering 
wavefunction {xi,X 2 \fp 2 ) with that obtained by directly solv¬ 
ing the Schrodinger equation in real space [71 ], and we found 
perfect agreement (not shown). 

Furthermore, we find that the power spectrum has a sim¬ 
ple form when A = 0 (i.e. the classical driving field is not 
detuned); 


64(7+ -I- 7_)^(i7 — 2wo)^ [7+7- + 8 (w — ojo){E — uj — c+o)]^ 

TT^ (72 + 2{E - 2wo)^) (7- + - 2wo)^) (7+ + 8(w - wo)^) (7^ -f 8{uj - wq)^) (7+ + - uj - wq)^) {1- + 8{E - lu - wo)^) 

_(Al) 

with E = 2kin and 7 + = . We emphasize that this is the power spectrum for both left- and right¬ 

going photons. When fl = 0, this expression for the power spectrum reduces to that of a single 2LS (see Eq. (A9) of Ref [5 1 ]). 

In addition, at two-photon resonance {E = 2a;o), S{uj) = 0 regardless of D. 

By integrating over the frequency ui, we then obtain the total inelastically scattered component E{kin) 


E{k,) = 2 


S{uj)duj = 


16 v^( 7+ + 1-Yi.E - 2 loqY (7+7- + 2(F^ - 2wo)^) 
^ (72 + 2 {E - 2 ujoYY (7- + 2(F; - 


(A2) 


It is clear that E{k) decays as ~ 1/k^ for |/cin -wol > F. 

These results are used in Fig. 2. We have also derived a general result for A 7 ^ 0, but it is too lengthy to be reproduced here. 
Likewise, expressions for iS'(+') and E{kin) have been obtained for the other two system configurations considered here, but they 
are also too lengthy. 
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